function [rate, azimuth] = SpeedAzimuth(ww, SHEI, FAI, longit, latit)  
	SHEI = SHEI*(pi/180);
	FAI = FAI*(pi/180);
	latit = latit*(pi/180);
	longit = longit*(pi/180);
	a = 6371;
    
	Vns = a*ww*cos(SHEI)*sin(longit-FAI);
	Vew = a*ww*(cos(latit)*sin(SHEI)-sin(latit)*cos(SHEI)*cos(longit-FAI));

	rate = ((Vns)^(2) + (Vew)^(2))^(1/2); 
	azimuth = 90 - (atan((Vns)/(Vew))*(180/pi));
end

%latitude 入纬度 经度longitude U 欧拉SHEN 欧拉FAI